Search for: All records

Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofR_B≈ 2 × 10⁵GM_•/c², whereM_•is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r⁻¹rather thanr^−3/2and the accretion rate$\dot{M}$is consequently suppressed by over 2 orders of magnitude below the Bondi rate${\dot{M}}_{\mathrm{B}}$. We find continuous energy feedback from the accretion flow to the external medium at a level of$\sim {10}^{-2}\dot{M}{c}^2\sim 5\times {10}^{-5}{\dot{M}}_{\mathrm{B}}{c}^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

 

The recent Chandra-JWST discovery of a quasar in thez≈ 10.1 galaxy UHZ1 reveals that accreting supermassive black holes were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity ofL_bol∼ 5 × 10⁴⁵erg s⁻¹, which corresponds to an estimated black hole (BH) mass of ∼4 × 10⁷M_⊙, assuming accretion at the Eddington rate. JWST NIRCAM and NIRSpec data yield a stellar mass estimate for UHZ1 comparable to its BH mass. These characteristics are in excellent agreement with prior theoretical predictions for a unique class of transient, high-redshift objects, overmassive black hole galaxies (OBGs) by Natarajan et al., that harbor a heavy initial black hole seed that likely formed from the direct collapse of the gas. Given the excellent agreement between the observed multiwavelength properties of UHZ1 and theoretical model template predictions, we suggest that UHZ1 is the first detected OBG candidate. Our assertion rests on multiple lines of concordant evidence between model predictions and the following observed properties of UHZ1: its X-ray detection and the estimated ratio of the X-ray flux to the IR flux, which is consistent with theoretical expectations for a heavy initial BH seed; its high measured redshift ofz≈ 10.1, as predicted for the transient OBG stage (9 <z< 12); the amplitude and shape of the detected JWST spectral energy distribution (SED) between 1 and 5μm, which is in very good agreement with simulated template SEDs for OBGs; and the extended JWST morphology of UHZ1, which is suggestive of a recent merge and is also expected for the formation of transient OBGs. As the first OBG candidate, UHZ1 provides compelling evidence for the formation of heavy initial seeds from direct collapse in the early Universe.

 

Supermassive black hole binary systems (SMBHBs) should be the most powerful sources of gravitational waves (GWs) in the universe. Once pulsar timing arrays (PTAs) detect the stochastic GW background from their cosmic merger history, searching for individually resolvable binaries will take on new importance. Since these individual SMBHBs are expected to be rare, here we explore how strong gravitational lensing can act as a tool for increasing their detection prospects by magnifying fainter sources and bringing them into view. Unlike for electromagnetic waves, when the geometric optics limit is nearly always valid, for GWs the wave-diffraction-interference effects can become important when the wavelength of the GWs is larger than the Schwarzchild radius of the lens, i.e.,$M_{\mathrm{lens}}\sim {10}^8{\left(\frac{f}{\mathrm{mHz}}\right)}^{-1}{M}_{\odot }$. For the GW frequency range explored in this work, the geometric optics limit holds. We investigate GW signals from SMBHBs that might be detectable with current and future PTAs under the assumption that quasars serve as bright beacons that signal a recent merger. Using the black hole mass function derived from quasars and a physically motivated magnification distribution, we expect to detect a few strongly lensed binary systems out toz≈ 2. Additionally, for a range of fixed magnifications 2 ≤μ≤ 100, strong lensing adds up to ∼30 more detectable binaries for PTAs. Finally, we investigate the possibility of observing both time-delayed electromagnetic signals and GW signals from these strongly lensed binary systems—that will provide us with unprecedented multi-messenger insights into their orbital evolution.

 

We look for simulated star-forming linear features such as the one recently discovered by van Dokkum et al. in the cosmological hydrodynamical simulationASTRID. Among the runaway black holes inASTRID, none are able to produce clear star-forming wakes. Meanwhile, flyby encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of black holes), reproduce remarkably well many of the key properties (length and linearity, recent star formation, etc.) of the observed star-forming linear feature. We predict that the feature will persist for approximately 100 Myr in such a system and hence constitute a rare event. The feature contains a partly stripped galaxy (withM_gal= 10⁹–10¹⁰M_⊙) and a dual black hole system (M_BH= 10⁵–10⁷M_⊙) in its brightest knot. The X-ray emission from AGN in the knot should be detectable in such systems. After 100–200 Myr from the first flyby, the galaxies merge, leaving behind a triple black hole system in a (still) actively star-forming early-type remnant of mass ∼5 × 10¹⁰M_⊙. Follow-up JWST observations may be key for revealing the nature of these linear features by potentially detecting the older stellar populations constituting the bright knot. Confirmation of such detections may therefore help discriminate a flyby encounter from a massive black hole wake to reveal the origin of such features.

 

Abstract Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ∼10–100 Hz, over a characteristic jet activity timescale ∼10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ∼10 events yr −1 . The detection rate drops to ∼1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events. 

While supermassive black-hole masses have been cataloged across cosmic time, only a few dozen of them have robust spin measurements. By extending and improving the existing Event Horizon Telescope (EHT) array, the next-generation Event Horizon Telescope (ngEHT) will enable multifrequency, polarimetric movies on event-horizon scales, which will place new constraints on the space-time and accretion flow. By combining this information, it is anticipated that the ngEHT may be able to measure tens of supermassive black-hole masses and spins. In this white paper, we discuss existing spin measurements and many proposed techniques with which the ngEHT could potentially measure spins of target supermassive black holes. Spins measured by the ngEHT would represent a completely new sample of sources that, unlike pre-existing samples, would not be biased towards objects with high accretion rates. Such a sample would provide new insights into the accretion, feedback, and cosmic assembly of supermassive black holes. 

We present a new method to simultaneously and self-consistently model the mass distribution of galaxy clusters that combines constraints from strong lensing features, X-ray emission, and galaxy kinematics measurements. We are able to successfully decompose clusters into their collisionless and collisional mass components thanks to the X-ray surface brightness, as well as use the dynamics of cluster members, to obtain more accurate masses exploiting the fundamental plane of elliptical galaxies. Knowledge from all observables is included through a consistent Bayesian approach in the likelihood or in physically motivated priors. We apply this method to the galaxy cluster Abell S1063 and produce a mass model that we publicly release with this paper. The resulting mass distribution presents different ellipticities for the intra-cluster gas and the other large-scale mass components as well as deviation from elliptical symmetry in the main halo. We assess the ability of our method to recover the masses of the different elements of the cluster using a mock cluster based on a simplified version of our Abell S1063 model. Thanks to the wealth of mutliwavelength information provided by the mass model and the detected X-ray emission, we also found evidence for an ongoing merger event with gas sloshing from a smaller infalling structure into the main cluster. In agreement with previous findings, the total mass, gas profile, and gas mass fraction are all consistent with small deviations from the hydrostatic equilibrium. This new mass model for Abell S1063 is publicly available, as the lenstool extension used to construct it.

 

We use the Galaxy Morphology Posterior Estimation Network (GaMPEN) to estimate morphological parameters and associated uncertainties for ∼8 million galaxies in the Hyper Suprime-Cam Wide survey withz≤ 0.75 andm≤ 23. GaMPEN is a machine-learning framework that estimates Bayesian posteriors for a galaxy’s bulge-to-total light ratio (L_B/L_T), effective radius (R_e), and flux (F). By first training on simulations of galaxies and then applying transfer learning using real data, we trained GaMPEN with <1% of our data set. This two-step process will be critical for applying machine-learning algorithms to future large imaging surveys, such as the Rubin-Legacy Survey of Space and Time, the Nancy Grace Roman Space Telescope, and Euclid. By comparing our results to those obtained using light profile fitting, we demonstrate that GaMPEN’s predicted posterior distributions are well calibrated (≲5% deviation) and accurate. This represents a significant improvement over light profile fitting algorithms, which underestimate uncertainties by as much as ∼60%. For an overlapping subsample, we also compare the derived morphological parameters with values in two external catalogs and find that the results agree within the limits of uncertainties predicted by GaMPEN. This step also permits us to define an empirical relationship between the Sérsic index andL_B/L_Tthat can be used to convert between these two parameters. The catalog presented here represents a significant improvement in size (∼10×), depth (∼4 mag), and uncertainty quantification over previous state-of-the-art bulge+disk decomposition catalogs. With this work, we also release GaMPEN’s source code and trained models, which can be adapted to other data sets.

 

We present estimates for the number of supermassive black holes (SMBHs) for which the next-generation Event Horizon Telescope (ngEHT) can identify the black hole “shadow”, along with estimates for how many black hole masses and spins the ngEHT can expect to constrain using measurements of horizon-resolved emission structure. Building on prior theoretical studies of SMBH accretion flows and analyses carried out by the Event Horizon Telescope (EHT) collaboration, we construct a simple geometric model for the polarized emission structure around a black hole, and we associate parameters of this model with the three physical quantities of interest. We generate a large number of realistic synthetic ngEHT datasets across different assumed source sizes and flux densities, and we estimate the precision with which our defined proxies for physical parameters could be measured from these datasets. Under April weather conditions and using an observing frequency of 230 GHz, we predict that a “Phase 1” ngEHT can potentially measure ∼50 black hole masses, ∼30 black hole spins, and ∼7 black hole shadows across the entire sky. 

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models, including a constant electron to magnetic pressure (constant βe model) with a population of non-thermal electrons as well as an R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations are enhanced with every pair added. On the contrary, we show that, at higher frequencies, a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of a large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function.